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This  thesis  presents  a  method  of  segmenting  unformatted 
alpha-numeric  characters.  Reconstructing  the  magnitude 
of  the  Fourier  transform  of  a  template  character  with  the 
phase  of  a  string  of  unformatted  characters  containing 
the  template  character  causes  all  characters  that  do  not 
have  the  magnitude  of  the  template  to  be  attenuated  in 
the  visual  domain.  The  template  character  will  not  be 
attenuated,  since  it  has  both  proper  magnitude  and  phase, 
and  a  peak  detector  can  find  the  most  probable  location 
of  the  character.  This  process  also  gives  a  first  choice 
for  what  the  character  is  (the  template).  Results  are 
presented  for  most  of  the  English  alphabet. 
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Many  U.  S.  Government  agencies  process  thousands  of  pages 
of  printed  material  each  month.  To  store  this  information 
in  easily  retrievable  and  processable  form,  it  has  to  be 
entered  into  a  computer.  Presently,  typists  perform  this 
task,  but  if  a  machine  could  read  from  these  sources  without 
human  intervention,  then  the  input  process  could  go  at 
machine  speeds,  which  are  much  faster.  This  could  potentially 
save  a  large  amount  of  human  labor. 

BaQ.kgr.Qund. 

All  current  Optical  Character  Recognition  (OCR)  techniques 
require  prior  knowledge  of  the  location  and  form  of  each  char¬ 
acter.  Many  useful  OCR  methods  are  available,  and  they  work 
reasonably  well  if  the  characters  are  typewritten  or  carefully 
handwritten  in  pre-printed  boxes  [1-4],  Most  existing  text 
is  not  so  formatted;  some  examples  are: 

a.  Books 

b.  Journals 

c.  Magazines 

d.  Other  published  material 

e.  ZIP  Codes  or  complete  addresses 

"Segmentation"  is  the  term  used  to  describe  the  process  of 
separating  the  characters  from  each  other  in  a  line  of  print. 
Most  people  have  seen  order  forms  that  can  be  machine-read  if 
characters  are  carefully  placed  in  pre-printed  boxes.  (See 
figure  1.1.)  This  is  one  form  of  segmentation.  But  this 


type  of  segmentation  is  not  available  when  machine-reading 
from  the  examples  listed  above.  An  OCR  technique  that  did 
not  require  the  input  characters  to  be  specifically  formatted 
for  OCR  purposes  would  be  much  more  useful. 

Problem. 

There  were  had  two  research  goals.  The  first  was  to 
develop  software  to  perform  Horev's  algorithm  [5]  on  the 
Data  General  Corporation  computers  in  the  Signal  Processing 
Laboratory.  This  goal  was  accomplished.  The  second  was  to 
test  the  feasibility  of  this  approach  for  segmenting  unfor¬ 
matted  alpha-numeric  characters.  This  report  contains  the 
results  of  the  feasibility  study. 


Scope . 

This  thesis  project  progressed  through  three  stages:  (1) 
background  research,  (2)  implementation,  and  (3)  analysis  of 
results.  The  background  research  was  limited  to  studying 
the  target-location  method  developed  by  Moshe  Horev.  His 
Master's  thesis  [5]  desribes  a  method  of  finding  a  target  in 
a  cluttered  scene.  Segmentation  of  characters  can  be  thought 
of  as  finding  a  target  (a  character  location)  in  a  cluttered 
scene  (the  line  of  print).  During  implementation,  the  study 
was  limited  to  English  alpha-numeric  characters,  since  results 
should  be  extendable  to  all  shapes  of  characters  and  numbers. 
Hand-drawn  characters  were  often  used  so  that  no  template 
would  exactly  match  the  character  in  the  input  scene,  and 
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Sample  OCR  Form  with  Segmentation  Boxes 
Figure  1.1 
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so  that  the  characters  could  be  made  different  sizes  to 


study  a  size  effect  that  is  discussed  in  the  next  chapter. 

The  results  are  presented  visually  in  this  report  to  clearly 
show  the  results  and  simplify  the  analysis. 

Assumptions . 

A  character  is  segmented  if  it  is  bounded  from  left  to 
right.  It  is  assumed  that  a  threshold  program  can  put  a 
"window"  around  the  segmented  character  if  its  left-most  and 
right-most  boundaries  are  not  removed  by  the  process.  Also, 
the  height  of  the  characters  are  assumed  to  be  set  to  a 
given  height  (determined  by  a  line-location  process)  before 
this  segmentation  algorithm  is  used  [2:1530], 

flyer_V-ie.w  si£  Chapters. 

The  next  chapter  contains  an  explanation  of  the  theory 
behind  Horev’s  algorithm,  with  emphasis  on  the  two-dimen¬ 
sional  discrete  Fourier  transform  ( 2D-DFT) .  The  third  chapter 
contains  a  general  description  of  the  software  development. 
Examples  and  an  interpretation  of  the  results  are  given  in 
chapter  four.  Finally,  chapter  five  gives  my  conclusion  and 
recommendations. 
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II.  Reconstruction  of  the  Fourier  Magnitude  of 
One  Scene  with  the  F our ier  Phase  of  Another 

General. 

To  process  any  visual  information,  a  domain  has  to  be 
chosen  that  supports  a  useful  calculation.  Some  researchers 
work  solely  in  the  spatial  domain — the  computations  are  per¬ 
formed  on  numerical  values  of  brightness  at  discrete  points 
across  the  two-dimensional  picture.  These  points  are  called 
pixels  (from  "picture  elements").  Other  researchers  [6,71 
have  suggested  that  the  human  visual  system  may  perform  its 
calculations  in  the  spatial  frequency  domain,  and  logically, 
the  greatest  amount  of  success  would  occur  if  machines  also 
processed  in  the  frequency  domain.  The  discrete  Fourier 
transform  (DFT)  is  used  to  compute  the  real  and  imaginary 
frequency  components  of  the  pixels.  AFIT  student  Moshe  Horev 
designed  an  algorithm  to  locate  targets  in  a  cluttered  scene, 
performing  his  calculations  in  the  frequency  domain  [5:9]. 
This  thesis  applies  a  much-simplified  version  of  his  method 
to  the  problem  of  character  segmentation.  Although  the  com¬ 
putation  occurs  in  the  frequency  domain,  the  results  will  be 
presented  in  the  visual  (spatial)  domain. 

In  two-dimensional  functions  (such  as  pictures),  the 
phase  of  the  signal  contains  most  of  the  form  information. 

The  importance  of  this  will  be  discussed  at  the  end  of  this 
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chapter;  first,  a  general  description  of  the  mathematics 
used  in  this  segmentation  algorithm  will  be  given. 


Two-Dimensional  Fourier  Transform. 

To  transform  the  discrete  pixels  from  the  visual  domain 
to  the  frequency  domain,  a  two-dimensional  discrete  Fourier  ‘ 
transform  (2D-DFT)  is  used.  The  pertinent  equation  for  the 
2D-DFT  is  [8:1173: 

M-1  N-1  2nxm  2nyn 

F(x,y)  =  V  V  f(m,n)  exp[-j( - )]exp[-j( - )]  2.1 

m=t)  n  =  0  M  N 

between  0  n  N-1  and  0  ^  m  ^  M-1.  Here,  f(m,n)  is  the 
brightness  value  of  the  pixel  at  location  (m,n)  where  m  and 
n  are  summed  over  the  area,  and  F(x,y)  is  the  Fourier  trans¬ 
form  at  any  location  (x,y). 

This  summation  can  be  simplified  if  the  dimensions  N  and 
M  are  equal  and  are  powers  of  two,  since  a  recursive  approach 
can  be  efficiently  used  to  perform  the  Fourier  transform. 
Specifically,  the  algorithm  I  used  [written  by  Ronald  W. 
Schaefer,  30  June  1978]  operates  only  on  a  square  array  with 
dimensions  which  are  powers  of  two.  It  performs  N  one-dimen¬ 
sional  Fourier  transformations  on  the  N  rows  of  the  two- 
dimensional  array,  then  transposes  the  array,  and  finishes 
by  again  performing  N  one-dimensional  Fourier  transforms. 

The  result  is  a  2D-DFT  [8:320-321]  that  is  transposed. 

Similarly,  the  inverse  two-dimensional  Fourier  transform 
(2D-IDFT)  is  computed  by  this  equation: 
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1  X-1  Y-1  2nmx  2nny 

f(m,n)  =  —  V  y  F{x,y)  exp[j( - )]exp[j( - )]  2.2 

XY  xTb  yrb  X  Y 

between  0  ^  x  ^  X-1  and  0  £  y  £  Y-1. 

Computation  of  Magnitude  and  Phase. 

The  magnitude  of  the  complex  Fourier  transform  was  cal¬ 
culated  by  the  standard  equation  (eqn.  2.3): 

2  2  v, 

MAG  =  [  REAL  +  IMAG  ]z  2.3 

i  i  i 


Similarly,  the  phase  was  calculated  by  equation  (2.4): 

-1  IMAGi 

ANG  =  tan  [ - ]  2.4 

i  REALi 

except  when  MAG,  =  0,  in  which  case  ANG.  =  0. 

i  l 

Reconstruction  of  the  Template  Magnitude  with  the  Phase 
&  ScQOg. 

From  Horev’s  work,  tho“key  to  segmentation  is  reconstruct¬ 
ing  the  magnitude  of  a  template  with  the  phase  of  a  scene; 
this  creates  a  new  complex  scene  where  segmentation  is 
possible  by  a  threshold  algorithm.  The  real  component  of 
the  new  complex  number  is  (equation  2.5): 

REAL  =  ENM  «  cos  (ANG  )  2.5 

i  i  i 


and  the  imaginary  component  is  (equation  2.6): 

IMAG  =  ENM  «  cos  (ANG  )  2.6 

i  i  i 
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where  ENM^  is  the  energy-normalized  value  of  MAG^.  The  norm¬ 
alization  factor  is  the  square  root  of  the  sum  of  the  squares 
of  each  magnitude  component.  By  repeated  application  of 
equations  2.5  and  2.6,  a  complex  file  is  recreated  from  a 
magnitude  file  and  a  phase  file.  Note  that  if  the  magnitude 
and  phase  of  the  same  scene  were  recombined  in  this  manner, 
the  original  (complex)  scene  would  be  recreated,  but  the 
pixels  would  be  energy-normalized. 

Importance  of  Phase  in  a  Scene. 

Research  has  shown  that  most  of  the  form  information  in  a 
scene  is  contained  in  the  phase  of  its  frequency  components 
[9].  For  a  typical  scene,  a  recognizeable  version  of  the 
scene  is  reproduced  if  a  file  of  unity  magnitudes  is  recon¬ 
structed  with  the  phase  of  the  scene  (and  then  inverse 
Fourier-transformed;  see  figure  2.1).  Having  the  whole  scene 
present  in  the  reconstructed  scene  will  not  aid  segmentation; 
the  reconstructed  scene  will  practically  match  the  original 
scene.  However,  the  size  of  a  character  relative  to  the 
2D-DFT  window  size  changes  the  amount  of  this  effect. 

Segmentation  can  be  accomplished  if  the  characters  in  the 
2D-DFT  window  are  relatively  large.  Normally,  the  most  energy 
in  a  frequency  spectrum  is  clustered  around  the  DC  term,  that 
is,  the  energy  is  in  the  low  frequency  components.  When 
unity  magnitudes  are  reconstructed  with  the  phase  of  a  scene, 
the  higher  frequencies  tend  to  be  emphasized.  This  will 
be  seen  visually  as  the  edges  in  a  scene.  When  a  larger 


character  is  in  the  window,  it  contains  mostly  low-frequency 
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Scene  Reconstructed  from  Ihase  with  Jnity  Magnitude 

Figure  2.1 
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components  (each  stroke  in  the  character  will  have  more  pixel 
to  represent  it  and  they  will  be  nearly  the  same  level  of 
black;  thus  there  will  be  few  high  frequency  changes).  The 
low-frequency  components  are  not  emphasized  when  recon¬ 
structed  with  unity  magnitudes,  and  this  scene  will  not  be 
recognizeable.  For  the  character  to  be  visible,  it  will  also 
need  its  magnitude  components.  Thus,  segmentation  is  more 
likely  to  be  possible  if  the  character  size-to-window  size 
ratio  is  large  enough  to  keep  the  phase  from  accurately 
representing  the  scene.  Furthermore,  in  the  actual  recon¬ 
struction  process,  the  magnitude  file  has  most  of  its  energy 
in  the  lower  frequency  components.  This  will  cause  the  high 
frequency  components  of  the  target  character  to  be  somewhat 
de-emphasized  back  to  their  normal  levels. 

In  this  research,  recognizeable  characters  were  still 
present  in  phase-only  scenes  when  one  of  the  characters 
already  filled  75%  of  the  window  (in  height;  i.  e.,  about  50% 
of  the  window  area).  To 'Fimpirify~ThTs  work,  only  large  char¬ 
acters  (total  character-to-window  ratio  above  75%)  were  used, 
except  when  working  with  original  characters  that  were  too 
small  to  optically  enlarge  to  that  ratio.  The  next  chapter 
describes  all  of  the  equipment  and  software  used  to  perform 
this  reconstruction  algorithm. 
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General. 

A  Signal  Processing  Laboratory  (SPL)  has  been  developed 
in  the  Electrical  Engineering  department  at  AFIT.  To  most 
effectively  use  this  resource,  software  must  be  developed  to 
process  digitized  signals  (audio  and  video,  although  this 
thesis  includes  only  video).  As  a  foundation  for  further 
video-processing  software,  this  chapter  contains  three 
sections:  the  first  section  is  a  general  description  of  the 
hardware  available  in  the  SPL  at  the  time  of  this  writing. 
The  second  describes  the  necessary  details  of  digitization 
and  processing  so  that  the  reader  has  a  grasp  of  what  is 
necessary  to  write  software  to  process  video.  Finally,  the 
third  section  will  briefly  describe  the  video-processing 
software  that  was  developed  concurrently  in  the  SPL.  Source 
listings  of  the  programs  are  given  in  the  Appendix. 

Hardware. 

The  first  step  in  understanding  how  video  information  is 
processed  is  to  gain  a  working  knowledge  of  the  equipment 
currently  available  in  the  SPL.  The  principal  components  of 
the  system  are  two  minicomputers  and  a  microcomputer.  They 
are: 


a.  Data  General  Nova  2  processor 

b.  Data  General  Eclipse  S/250  processor 

c.  Cromemco  Z-2D  microcomputer 
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The  two  Data  General  machines  are  used  for  the  actual  data 
processing,  while  the  Cromemco  computer  is  a  Z-80  based 
processor  that  controls  the  analog- to-digital  and  digital- 
to-analog  conversion  devices.  The  A/D  and  D/A  boards  were 
built  by  Tecmar,  Inc.  and  they  have  the  capability  of  digit¬ 
izing  pictures  of  64 ,  128,  256,  or  512  pixels  per  line  by  64, 
128,  or  256  lines.  Input  is  from  a  standard  studio  vidicon 
camera  (Cohu  Electronics,  Inc.  model  6150-000),  and  output 
is  to  any  standard  monitor  (or  television  set,  if  the  signal 
is  modulated  onto  a  standard  carrier  frequency). 

The  Nova  computer  starts  and  stops  the  Cromemco  via  a 
fortran  subroutine  call  named  CHANNEL  [10].  This  subroutine 
performs  two  tasks;  one  is  to  pass  proper  parameters  to  the 
Cromemco  for  signal  handling,  and  the  other  is  to  retrieve  or 
send  the  actual  data.  The  Z-2D  is  an  eight  bit  processor, 
while  both  Data  General  machines  have  sixteen  bit  words. 

The  subroutine  call  performs  the  necessary  conversion  to 
interface  the  machines  and  writes  or  reads  directly  to  the 
disk  drives.  (The  SPL  currently  has  a  ten  megabyte  and  two 
twenty  megabyte  Data  General  disk  drives.)  Once  the  digit¬ 
ized  picture  is  stored  on  disk,  either  the  Nova  or  the 
Eclipse  can  access  the  information  and  process  it.  However, 
the  Eclipse  cannot  directly  communicate  with  the  Cromemco; 
it  can  only  access  the  disk  drives  on  which  the  Nova  stores 
the  video  data. 

The  steps  necessary  to  handle  video  signal  processing  are: 


Input:  a.  Activate  vidicon,  Cromemco,  and  Nova 
systems 

b.  Call  subroutine  CHANNEL  from  Nova  to 
operate  A/D  and  transfer  data  from 
Cromemco  to  Nova 

c.  Write  video  data  to  disk 

Process  data:  a.  Nova 

b.  Eclipse 

Output:  a.  Insure  that  video  data  is  on  a  disk 
that  is  accesible  to  the  Nova 

b.  Activate  Cromemco,  Nova,  and  output 
device  (e.  g.  monitor) 

c.  Call  subroutine  CHANNEL  to  transfer 
data  from  Nova  to  Cromemco  for  D/A 

Note  that  the  Cromemco  computer  is  completely  transparent  to 

the  Nova;  the  user  simply  calls  the  subroutine  and  the  Nova 

automatically  either  inputs  the  video  frame  that  the  Cromemco 

has  ready  or  outputs  a  video  frame  to  the  Cromemco. 


Video  Digitization  and  Processing . 

The  software  running  in  the  Cromemco  and  switches  on  the 
video  digitization  boards  set  certain  constraints  on  the  form 
of  the  video  data.  The  picture  frame  is  usually  digitized 
at  256  by  256  pixels.  Furthermore,  the  Tecmar  A/D  board 
digitizes  the  picture  into  sixteen  gray  levels  (zero  corre¬ 
sponds  to  black,  and  fifteen  corresponds  to  white).  Thus, 
four  bits  are  needed  to  represent  one  pixel.  The  A/D  con¬ 
verter  outputs  two  pixels  at  a  time,  as  an  eight  bit  word. 
Since  the  Data  General  computers  use  sixteen  bit  words,  the 
fortran  subroutine  that  activates  the  Cromemco  reads  two  of 
the  eight  bit  words  and  creates  a  sixteen  bit  word  that  is 
written  to  disk.  What  results  is  a  sixteen  bit  word  with 
four  pixels  stored  in  it.  (See  figure  3.1.) 


If  any  operation  is  to  be  performed  on  the  pixel's  sep¬ 
arately,  they  first  have  to  be  "unpacked"  from  the  sixteen 
bit  word.  After  any  processing,  the  four  adjacent  pixels 
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must  be  "repacked"  if  the  video  scene  is  to  be  displayed  on 
a  standard  output  device.  Also  note  that  the  packed  version 
of  a  video  file  is  an  efficient  way  to  store  the  picture. 

Even  though  the  frame  is  best  stored  as  packed  pixels, 
after  processing,  the  pixels  are  often  no  longer  four  bit 
integers,  ror  example,  after  Fourier  transformation,  the 
pixels  are  complex  numbers;  their  magnitudes  are  real  numbers 
often  greater  than  15.  Thus,  in  general,  it  is  easier  to 
store  transformed  versions^  of  a  scene  in  an  unpacked  form. 

The  algorithm  presented  in  this  thesis  is  actually  seven 
stand-alone  programs.  There  are  three  major  reasons  for  this 
First,  the  2D-DFT  and  2D-IDFT  programs  were  supplied  as  sep¬ 
arate,  executable  programs.  This  immediately  separated  the 
algorithm  into  three  parts:  (1)  before  the  2D-DFT,  (2) 
between  the  2D-DFT  and  the  2D-IDFT,  and  (3)  after  the  2D-IDFT 
Second,  core  size  limitations  restricts  the  length  of  many 
programs.  Finally,  many  steps  of  the  algorithm  have  to  be 
executed  only  once.  For  example,  once  the  phase  of  a 


scene  is  stored  (or  the  magnitude  of  a  template),  those  data 
do  not  need  to  be  created  again.  Therefore,  the  algorithm 
is  best  executed  by  separate  programs. 

Software  Development. 

Figure  3.2  contains  a  flow  chart  of  the  programs  that  make 
up  this  segmentation  algorithm.  The  development  of  each  pro¬ 
gram  (except  DIRECT  and  INVERSE,  which  were  supplied  programs) 
is  described  in  general  terms  in  the  next  paragraphs. 

The  first  of  the  eight  programs  is  called  "VIDEO."  This 
program  is  a  general-purpose  video  I/O  program  that  sets  up 
the  parameters  for  CHANNEL — the  subroutine  that  activates 
the  Cromemco  computer.  Certain  parameters  are  always  con¬ 
stant  for  video  I/O,  but  others  must  be  given  by  the  user: 
input  or  output,  monitor  display  time,  and  disk  filename. 
Program  VIDEO  requests  these  and  then  calls  subroutine 
CHANNEL. 

The  first  program  that  actually  performs  any  processing 
on  the  data  file  is  called  "VTOC"  (Video  TO  Complex).  This 
program  performs  three  steps:  (1)  It  reads  from  a  video  file 
and  unpacks  the  pixels;  (2)  It  "blanks  out"  a  portion  of  the 
picture  if  the  user  requests  this  operation;  and  (3)  It 
writes  the  pixels  out  to  disk  as  complex  numbers  (imaginary 
part  zero).  The  blanking  step  is  necessary  since  the  video 
digitizer/Cromemco  combination  does  not  completely  digitize 
a  256  by  256  picture.  Twelve  pixels  on  the  left  end  of  each 
line  and  the  last  23  lines  are  incorrectly  digitized.  See 
figure  3«3.  Since  this  noise  will  add  unwanted  frequency 
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components  to  the  Fourier  transform,  the  noise  was  blanked 
to  a  constant  level.  The  video-to-complex  transformation  is 
needed  so  that  the  video  file  is  readable  by  the  2D-DFT  pro¬ 
gram.  This  program  (DIRECT)  reads  and  writes  complex  numbers 
only.  Therefore,  the  only  preprocessing  necessary  is  the 
transformation  from  video  to  complex  numbers.  (Plus  blank-' 
ing,  if  the  user  desires.) 

After  Fourier  transformation,  the  magnitude  of  the  template 
must  be  calculated,  as  must  the  phase  of  the  scene.  A  pro¬ 
gram  to  perform  these  operations  is  "OREAD"  (for  Complex 
READ).  This  program  computes  the  magnitude  (see  equation 
2.3)  or  phase  (see  equation  2.4)  of  a  complex  file  and  writes 
the  output  as  real  numbers  to  a  new  disk  file.  At  this  point, 
the  complex  files  need  not  be  saved  — no  more  operations  are 
performed  on  them,  and  the  original  video  file  is  the  most 
easily  understood  representation  of  the  scene.  Their  disposal 
results  in  a  large  savings  of  disk  space. 

The  next  program  performs  a  simplified  version  of  Horev's 
algorithm;  the  program  is  named  "HRV."  It  performs  the 
reconstruction  of  a  magnitude  file -with  a  phase  file, 
creating  a  complex  file  that  the  2D-IDFT  (INVERSE)  can  read. 
Equations  2.5  and  2.6  are  used  to  create  the  complex  terms. 

The  steps  performed  by  program  HRV  are:  (1)  Computation  of 
the  square  root  of  the  sura  of  the  squares  of  each  magnitude 
component  (for  energy  normalization),  (2)  Reconstruct  the 
new  complex  scene  from  the  energy-normalized  magnitudes,  and 
(3)  Writing  the  new  complex  components  to  disk.  After  the 


reconstruction  operation  and  the  2D-IDFT,  the  data  file  is 
complex.  The  magnitude  of  this  complex  scene  contains  the 
visual  information,  so  the  scene  can  be  stored  much  more 
efficiently  in  magnitude  form.  Program  CREAD  calculates  the 
magnitude  of  this  file;  no  other  computation  is  performed  on 
the  complex  numbers. 

However,  since  segmentation  was  verified  in  this  work  by 
viewing  the  reconstructed  scene,  at  least  one  more  program 
had  to  be  written.  This  program  converts  the  (reconstructed) 
magnitude  file  to  repacked  video  pixels  for  display  on  a 
monitor.  Two  versions  of  this  type  of  program  were  written. 
The  first,  named  PIX1 ,  clips  the  magnitudes  into  a  zero  to 
fifteen  range.  This  is  necessary  since  the  magnitudes  of 
the  reconstructed  file  are  typically  of  the  order  0.0002  to 
2.0;  they  had  to  be  linearly  scaled  to  the  range  zero  to  fif¬ 
teen.  Equation  3.1  was  used  to  clip  the  magnitudes  to  the 
proper  range. 

PIXELi  =  15  *  C  ( Rmax  -  MAG  )/<Rmax  -  Rmin)]  3.1 

In  equation  3.1,  MAG^  is  the  magnitude  of  the  pixel  after 
the  CREAD,  Rmax  is  the  largest  magnitude,  and  Rmin  is  the 
smallest  magnitude.  Clipping  occurs  when  the  magnitude 
is  outside  of  a  user-specified  range.  If  the  magnitude  is 
above  the  range,  the  pixel  is  set  to  15.  If  the  magnitude 
is  below  the  range,  the  pixel  is  set  to  zero.  Clipping  the 
magnitudes  does  not  properly  perform  segmentation,  however, 
since  the  magnitudes  that  are  clipped  to  zero  will  be 


indistinguishable  from  the  black  (zero)  pixels  that  make 
up  the  segmented  character.  Thus,  PIX1  had  to  be  changed 
from  "clipping”  to  "exclusion";  pixels  outside  the  range 
were  not  printed  by  the  new  program.  Only  the  pixels  that 
fall  into  the  specified  range  are  printed.  If  the  proper 
range  has  been  selected,  then  only  the  segmented  character 
will  be  visible.  This  program  was  named  "PICK",  and  the 
results  presented  in  the  next  chapter  were  created  with 
program  PICK,  except  that  11  was  used  instead  of  15  in 
equation  3.1  so  that  there  would  be  a  pixel  level  "jump" 
from  the  white  background  (excluded  pixels)  to  the  character 
pixels. 

Another  useful  program  that  is  also  not  in  the  algorithm 
is  called  DISPLAY.  This  program  was  written  to  print  a  video 
scene  on  a  Printronix  P-300  lineprinter.  This  program  con¬ 
verts  the  pixels  to  arrays  of  lineprinter  dots  that  represent 
the  sixteen  gray  levels.  The  video  figures  included  in  this 
thesis- were  printed  from-^video  files  via  DISPLAY. 

With  this  software  developed,  different  sets  of  characters 
were  processed.  The  results  are  given  in  the  next  chapter. 
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IV.  Results  of  the  Reconstruction  Algorithm 


General. 

The  last  two  chapters  have  developed  the  theory  and  given 
a  brief  description  of  the  software  for  the  reconstruction 
algorithm.  This  chapter  contains  some  of  the  results 
achieved  when  segmenting  English  characters.  The  complete 
alphabet  is  not  used,  since  there  was  not  enough  time  to 
process  all  of  the  possible  two  or  three  letter  pairs.  Most 
of  the  characters  presented  were  handwritten  by  one  of  the 
students  in  the  SPL.  The  only  requirements  placed  on  the 
size  or  shape  of  the  characters  is  that  they  be  of  a  speci¬ 
fied  height  and  that  they  be  nonscript.  The  specified 
height  served  two  purposes:  first,  it  forced  the  characters 
to  fill  up  a  large  part  of  the  window  (as  described  on  pages 
8  and  10  in  chapter  2),  and  second,  it  was  assumed  that  an 
optical  zoom  feature  can_rbe  implemented  to  insure  that  all 
characters  are  of  approximately  the  same  height.  The  non¬ 
script  requirement  was  to  keep  the  style  of  each  character 
to  the  general,  standard  shape. 

Segmentation  Results. 

The  first  segmentation  results  presented  are  of  the  scene 
WKZ”  (figure  4.1).  The  phase  of  this  scene  was  calculated, 
and  then  the  magnitude  of  "Z"  (figure  4.2)  was  calculated. 
These  were  recombined,  inverse  Fourier  transformed,  and 
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then  the  magnitude  of  the  recombined  scene  was  saved.  Pro¬ 
gram  PICK  was  used  to  print  three  lower  ranges  of  pixel  levels. 
KZZHRV.V4  (figure  4.3)  has  a  maximum  of  0.2910  and  a  minimum 
of  0.0.  Recall  that  the  magnitude  range  after  the  inverse 
Fourier  transformation  is  0.0002  to  2.0.  (A  note  on  file 
notation:  the  ”.V4"  means  that  the  file  is  video  file  number 
four;  the  "HRV"  means  that  the  recombination  process  has 
been  completed;  the  character  in  front  of  "HRV"  is  the  tem¬ 
plate;  and  the  characters  in  front  of  the  template  character 
are  the  scene  characters.)  In  figure  4.3,  a  portion  of  the 
"K"  is  present,  while  the  "ZM  still  has  its  whole  width  pres¬ 
ent.  Dropping  to  a  lower  pixel  range  (figure  4.4)  removes 
more  of  the  "K” ,  while  the  MZW  width  is  still  completely 
present.  (The  maximum  magnitude  printed  is  0.2190.)  Dropping 
to  a  still  lower  maximum  magnitude,  KZZHRV.V2,  we  see  that 

the  "K"  is  virtually  removed,  while  the  left-to-right  range 

of  the  "ZM  is  easily  seen  (figure  4.5).  The  maximum  level  is 
0.1460,  with  the  minimum  still  0.0.  — 

Next,  "K"  was  used  as  a  template  to  remove  the  "Z"  (see 

figure  4.6).  KZKHRV.V3  (figure  4.7)  has  a  maximum  magnitude 
of  0.2120  and  a  minimum  of  0.0.  The  MZ f  is  virtually  gone. 
Dropping  to  a  maximum  of  0.1410  (figure  4.8)  or  0.07060  (fig¬ 
ure  4.9)  shows  that  the  "K"  is  still  present,  while  there  is 
no  width  to  the  WZP  at  all.  Even  if  the  remaining  MZ”  pixels 
of  KZKHRV.V1  (figure  4.9)  are  used  in  bounding  "K",  there  will 
be  little  error  in  the  placement  of  the  window  around  "K." 

The  characters  "K"  and  "Z**  both  have  a  lot  of  diagonal 


Reconstructed  Scene  KZZHRV.V4 
Figure  4.3 


Reconstructed  Scene  KZZHRV.V3 
Figure  4.4 
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Reconstructed  Scene  KZKHRV.V1 
Figure  4.9 
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information  in  them,  but  the  reconstruction  process  displays 
only  the  diagonal  information  that  is  contained  in  the  temp¬ 
late;  the  diagonal  in  "Zf  is  at  a  different  angle  than  that  in 
"K."  The  "Z"  is  segmented  mostly  by  its  horizontal  informa¬ 
tion,  while  the  "K"  is  segmented  by  its  vertical  information. 

The  next  set  of  reconstruction  results  is  with  the  letter 
pair  "MT"  (figure  4.10).  The  magnitude  of  the  template 
character  "M"  (figure  4.11)  was  calculated  and  used  to  segment 
the  "M"  from  the  "T;‘;  MTMHRV.V5  (figure  4.12)  has  a  maximum 
pixel  level  of  0.3370  and  a  minimum  pixel  magnitude  of  0.0. 

The  few  sparse  pixels  of  the  "T"  could  give  a  threshold 
algorithm  some  difficulty,  but  dropping  the  maximum  magnitude 
level  via  PICK  down  to  0.2700  gives  us  MTMHRV.V4  (figure  4.13). 
The  "T"  is  completely  removed.  To  be  completely  general,  "T" 
was  then  used  as  a  template  (figure  4.14)  to  remove  the  "M." 
Both  "M"  and  "T"  have  a  lot  of  vertical  information,  and 
MTTHRV.V5  (figure  4.15)  shows  that  all  of  the  vertical 
strokes  in  the  characters  are  present.  The  maximum  magnitude 
included  is  0.3440,  with  the  minimum  still  0.0.  Dropping 
the  maximum  level  to  0.2750  creates  MTTHRV.V4  (figure  4.16). 

The  crossbar  on  the  "T"  is  still  completely  present,  but 
the  right-hand  vertical  stroke  in  the  "M"  would  probably  be 
sensed  by  a  threshold  algorithm.  Dropping  the  maximum  to 
0.2070  in  MTTHRV.V3  (figure  4.17)  begins  to  remove  some  of 
the  crossbar  of  the  "T",  while  the  right-hand  bar  in  the  "M" 
is  still  present.  However,  in  either  figure  4.16  or  4.17, 
we  can  see  that  even  if  the  right-hand  vertical  pixels  of  the 
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Reconstructed  Scene  MTMHRV.V5 
Figure  4.12 
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Reconstructed  Scene  MTTHRV.V5 
Figure  4.15 


Reconstructed  Scene  MTTHRV.V4 
Figure  4.16 
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"M"  were  included  as  part  of  the  nT"  (by  thresholding),  the 
character  "T"  will  still  be  adequately  segmented. 

A  pair  of  characters  that  is  much  more  difficult  to  seg¬ 
ment  is  shown  in  figure  4.18.  Both  the  "W"  and  the  "A”  have 
many  similar  features.  To  further  stress  the  algorithm, 
serifs  were  included  on  the  template  "W"  (figure  4.19).  The 
peak  of  the  "A”  and  the  center  of  the  "W"  both  are  the  same: 
a  WAWHRV.V2  (figure  4.20)  has  a  maximum  level  of 

0.3670  (minimum  0.0),  and  the  "W  is  very  much  present;  the 
"A"  is  not  removed  enough  to  segment  the  "W" ,  however.  In 
the  next  figure  (4.21),  WAWHRV.V3  has  its  maximum  down  to 
0.3058,  and  the  "A"  is  now  just  two  points.  The  pixel  . 
levels  remaining  in  the  "W"  are  higher  than  those  in  the  "A", 
but  both  characters  may  still  meet  a  threshold. 

Using  an  WA"  as  a  template  (figure  4.22)  should  give  sim¬ 
ilar  results,  and  it  does.  In  figure  4.23  WAAHRV.V2  (max¬ 
imum  is  0.1480,  minimum  0.0)  shows  that  the  ”/n  is  located 
in  the  "A"  and  the  "V  is  found  in  the  'LWI* .  Dropping  to  a 
lower  maximum  magnitude  removes  much  of  the  "W"  (figure  4.24), 
but  both  characters  would  most  likely  meet  a  threshold 
criterion.  (Dropping  to  a  level  that  does  not  include  the 
"W"  also  does  not  include  the  "A”,  since  they  have  the  same 
maximum  magnitudes  before  conversion  to  pixels.)  The  maximum 
magnitude  was  dropped  to  0.0738.  If  both  characters  were 
included  in  the  segmentation  window,  the  correlation  process 
would  pick  out  the  "A." 

An  even  more  similar  pair  of  characters  is  ”0W  and  nQ.M 


Reconstructed  Scene  WAWHRV.V2 
Figure  4.20 


Reconstructed  Scene  WAWHRV.V3 
Figure  4.21 
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Reconstructed  Scene  WAAHRV.V1 
Figure  4.24 
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The  only  difference  between  these  two  characters  is  the  extra 


stroke  in  the  "Q."  Using  "0"  as  a  template  should  segment 
both  the  "0's"  and  "Q’s";  this  does  happen.  The  correlation 
process  that  occurs  after  segmentation  would  have  to  choose 
the  proper  character.  Figure  4.25  is  of  the  character  pair 
"QR."  Both  have  "roundness"  to  their  shapes  (the  top  half  of 
the  nRw  is  reasonably  round).  Using  "0"  as  the  template 
(figure  4.26),  the  ”Q ”  was  segmented  (figure  4.27).  QR0HRV.V2 
has  a  maximum  level  of  0.5100,  with  a  minimum  of  0.0.  A 
threshold  could  easily  pick  off  the  "Q"  and  not  the  "R." 
Dropping  down  to  a  lower  range  (figure  4.28),  QR0HRV.V1  has 
most  of  the  "R"  removed.  The  maximum  is  0.4420,  and  the 
minimum  is  0.0348.  Reducing  the  maximum  to  a  level  similar 
to  the  levels  used  in  the  examples  above  (maximum:  0.3741) 
results  in  QR0HRV.V3  (figure  4.29).  There  are  no  "R"  pixels, 
while  the  "Q"  is  still  present  on  the  left  and  right. 

Similarly,  many  characters  have  vertical  strokes  (like  "I") 
and  could  be  highlighted  by  using  "I"  as— a  template.  Other 
characters  have  two  vertical  strokes  (H,  M,  N,  and  U)  and 
using  "I"  as  a  template  should  pick  up  the  vertical  lines 
in  these  characters.  To  test  this  theory,  two  I's  were  used 
as  a  template  (figure  4.30)  to  segment  the  "U"  from  the  "S" 
in  the  scene  "US”  (figure  4.31).  USIIHRV. V 1  has  a  maxi¬ 
mum  magnitude  of  0.5973  and  a  minimum  of  0.0  (figure  4.32). 
Most  of  the  "S"  is  gone,  while  the  "U"  is  clearly  visible. 

A  threshold  should  be  able  to  pick  off  the  higher  level  of 
pixels  that  are  present  in  the  "U." 


However,  dropping  to  a 
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Reconstructed  Scene  QR0HRV.V2 
Figure  4.27 


Reconstructed  Scene  QR0HRV.V3 
Fifure  4.29 
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Template  II 
Figure  4.30 


Scene  US 
Figure  4.31 
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lower  maximum  magnitude  level  in  USIIHRV.V3  (figure  4.33), 
with  the  maximum  level  at  0.5347,  we  see  that  the  level  of 
the  "U"  is  far  above  the  level  of  the  "S."  From  these  last 
two  examples,  we  can  see  that  care  must  be  used  in  selecting 
the  template  characters. 

The  next  two  sets  of  characters  presented  both  have  the 
characters  touching  each  other  or  overlapping  vertically. 
This  is  where  segmentation  is  most  difficult  with  methods 
that  attempt  to  segment  characters  by  finding  gaps  between 
them;  these  methods  cannot  find  gaps  between  touching 
characters. 

The  first  set  is  the  three-character  group  "FAD"  (figure 
4.34).  The  "F"  touches  the  "A"  and  the  "A"  touches  the  "D." 
The  template  character  is  "F"  (figure  4.35).  FADFHRV . V3  has 
a  maximum  of  0.2710  and  a  minimum  of  0.0  (figure  4.36).  The 
"F”  can  easily  be  segmented  by  thresholding.  Dropping  the 
maximum  level  to  0.1450  in  FADFHRV. V2  (figure  4.37),  the  "F" 
is  still  quite  visible,  and  its -Width  would  be  easily  picked 
up  by  thresholding. 

The  second  set  of  harder-to-segment  characters  is  the 
three-letter  group  "GYJ"  (figure  4.38).  The  "G"  touches  the 
"Y"  and  the  "Y"  and  "J"  overlap  vertically.  The  template 
character  is  "Y"  (figure  4.39).  GYJYHRV.V6  is  shown  in 
figure  4.40.  Its  maximum  is  0.4070,  and  its  minimum  is 
0.0679.  The  "Y"  clearly  has  darker  pixel  levels  than  the 
"G";  dropping  the  maximum  to  0.3400  and  the  minimum  to 


0.0  results  in  GYJYHRV. V5  (figure  4.41).  The  "G"  is  virtually 
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Reconstructed  Scene  GYJYHRV.V5 
Figure  4.41 


removed,  and  the  "Y"  can  be  picked  up  by  a  threshold. 

In  the  examples  included  in  this  chapter,  large,  hand¬ 
written  characters  were  used.  With  this  algorithm,  smaller 
characters  were  not  segmentable;  when  smaller  characters 
were  used,  the  original  scene  was  recreated  in  virtually 
every  attempt.  For  example,  the  typewritten  scene  "The  quick 
(figure  4.42)  was  used  to  test  the  algorithm;  the  scene  has 
a  low  character  size-to-window  size  ratio.  The  template  was 
"T"  (figure  4.43).  The  segmentation  scene,  THE8THRV. V9  (fig¬ 
ure  4.44),  shows  that  most  of  the  characters  are  still  pres¬ 
ent;  a  window  formed  by  thresholding  would  not  reduce  the 
amount  of  characters  in  the  correlation  process.  This 
occurred  because  the  phase  information  contained  the  whole 
scene.  With  the  larger  characters,  however,  excellent 
results  occurred.  Only  with  the  set  "WA”  was  there  any 
problem.  This  shows  why  a  two-step  process  is  necessary  for 
unformatted  character  recognition:  the  segmentation  process 
will  only  bound  the  most  -likely  'location  of  a  character.  A 
correlation  process  will  be  needed  to  insure  that  the  char¬ 
acter  is  properly  recognized. 


Scene  "The  quick" 
Figure  4.42 
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Figure  4.43 
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V.  Conclusions  and  Recommendations 


.Re-View. 

Chapter  two  developed  the  theory  behind  this  segmentation 
algorithm.  Chapter  three  described  the  software  that  was 
developed  to  perform  the  reconstruction  algorithm,  and 
detailed  listings  are  included  in  Appendix.  Chapter  four 
presented  the  segmentation  results  for  half  of  the  English 
uppercase  alphabet.  Two  conclusions  can  be  reached,  and 
I  have  three  recommendations. 

Conclusions. 

A  simplified  version  of  Horev’s  algorithm  shows  great 
promise  in  segmenting  characters.  Since  characters  are 
usually  not  rotated,  and  their  height  can  be  set  to  a  con¬ 
stant  level  by  an  optical  process  before  segmentation,  the 
segmentation  algorithm  only  has  to  find  the  lef t-to-right 
range  of  each  character.  Although  researchers  [11]  have  shown 
that  Horev's  algorithm  has  certain  difficulties  performing 
its  task  (locating  targets  in  a  cluttered  scene  indepen¬ 
dent  of  size,  location,  or  rotation),  it  performs  well  in 
the  special  case  of  characters  on  a  contrasting  background. 

Regardless  of  this  algorithm’s  performance  as  shown  in 
chapter  four,  there  is  a  limitation.  The  algorithm  does  not 
perform  well  on  "small"  characters  (in  comparison  to  the 
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2D-DFT  window  size).  Either  the  characters  have  to  be  made 
larger  (optically  or  by  software),  or  the  window  has  to  be 
made  smaller  (a  software  change).  Either  or  both  of  these 
changes  could  conceivably  make  this  algorithm  work  for  most 
character  sizes. 

Recommendations . 

Three  recommendations  can  be  made.  They  are:  (1)  Software 
has  to  be  written  to  reduce  the  window  size  for  the  2D-DFT; 

(2)  An  algorithm  has  to  be  written  to  perform  the  threshold 
operation  and  place  a  "boundary”  around  the  segmented  char¬ 
acters);  and  (3)  The  effect  of  noise  on  the  algorithm  has 
to  be  studied. 

A  program  can  be  written  to  place  a  window  around  any 
location  in  the  scene.  A  256  by  256  window  is  used  in  this 
algorithm;  a  64  by  64  window  should  work  well  for  very  small 
characters  (e.  g.  characters  from  printed  material).  This 
will  also  reduce  the  2D-DFT  runtime  by  a  large  factor. 
Segmentation  of  smaller  characters  can  be  studied  if  this 
type  of  program  is  written.  The  main  reason  for  a  software 
change,  and  not  changing  the  video  digitizer  to  a  smaller 
size,  is  speed.  The  hardware  can  digitize  256  by  256  pixels 
practically  as  fast  as  64  by  64  pixels;  if  the  window  pro¬ 
gram  is  developed,  it  can  partition  the  256  by  256  pixel 
scene  into  16  subscenes  very  rapidly,  instead  of  requiring 
16  separate  input  operations. 
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An  algorithm  to  place  a  boundary  around  the  segmented 
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character  has  to  be  able  to  include  the  proper  pixel  magnitude 
range  that  represents  the  character  that  should  be  segmented. 
More  research  is  needed  to  develop  this  general  threshold 
criterion.  Table  I  lists  the  magnitudes  included  as  pixels 
in  the  scenes  in  this  thesis.  Since  all  the  parameters 
(except  what  each  character  is)  are  held  constant  (height  of 
the  characters,  lighting,  and  input  device),  the  segmented 
character  will  always  be  in  the  lowest  one-fifth  of  the  range 
of  magnitudes  in  the  reconstructed  scene  (See  Table  I).  Also, 
a  minimum  number  of  magnitudes  (or  pixels)  are  always  needed 
to  make  up  even  the  smallest  character.  A  threshold  algorithm 
would  have  to  consider  these  characteristics,  and  possibly 
more,  if  it  is  to  accurately  segment  characters  without  human 
intervention.  Once  the  algorithm  is  developed,  the  segmented 
characters  are  ready  to  be  read  by  a  correlation  process  [12], 
and  there  is  already  a  first  choice  for  each  character:  the 
template. 

The-  noise  was  not  removed  from  any  scene  (except  for  the 
noise  blanked  with  the  progam  VTOC;  see  chapter  3>  page  14) 
during  this  algorithm.  It  was  left  in  to  check  its  effect  on 
the  algorithm.  There  are  primarily  two  types  of  noise  that 
would  be  included  in  a  video  scene.  The  first  kind  of  noise 
is  from  the  digitization  process.  For  example,  if  the  char¬ 
acters  are  printed  on  certain  colors  of  paper,  there  may  not 
be  a  clear  distinction  between  the  background  and  the  char¬ 
acters.  (The  scenes  in  this  thesis  are  black  on  white,  and 
the  digitizer  rarely  digitized  the  white  background  as 
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Table  I 

Maximum  and  Minimum  Magnitude  Ranges  Included 
in  each  Scene  by  Program  PICK 


:gmentation 

Inclusion 

Magnitude 

Scene 

Maximum 

Minimum 

KZZHRV.V4 

0.2910 

0.0000 

KZZHRV.V3 

0.2190 

0.0000 

KZZHRV.V2 

0. 1460 

0.0000 

KZKHRV.V3 

0.2120 

0.0000 

KZKHRV.V2 

0.1410 

0.0000 

KZKHRV. V 1 

0.0706 

0.0000 

MTMHRV.V5 

0.3370 

0.0000 

MTMHRV.V4 

0.2700 

0.0000 

MTTHRV.V5 

0.3440 

0.0000 

MTTHRV.V4 

0.2750 

0.0000 

MTTHRV.V3 

0.2070 

0.0000 

WAWHRV.V2 

0.3670 

0.0000 

WAWHRV.V3 

0.3058 

0.0000 

WAAHRV.V2 

0.1480 

0.0000 

WAAHRV. V 1 

0.0738 

0.0000 

QR0HRV.V2 

0.5100 

0.0000 

QR0HRV.V1 

0.4420 

0.0348 

QR0HRV.V3 

—  __  rfeWI 

0.0000 

USIIHRV. V 1 

0.5973 

0.0000 

USIIHRV.V3 

0.5347 

0.0000 

FADFHRV.V3 

0.2710 

0.0000 

FADFHRV.V2 

0.1450 

0.0000 

GYJYHRV. V6 

0.4070 

0.0679 

GYJYHRV.V5 

0.3400 

0.0000 
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4  absolutely  white  pixels.)  The  second  noise  is  the  "frozen- 

frame"  effect  [4],  Random  noise  is  present  in  every  scene  that 
we  view,  but  it  is  not  seen  since  the  noise  is  not  correlated 
from  scene  to  scene.  But  when  the  scene  is  "frozen"  (digit¬ 
ized),  the  random  noise  is  clearly  visible. 

The  first  type  of  noise  can  often  be  removed  by  a  simple 
threshold  device  that  sets  all  pixels  below  a  certain  level 
to  black,  and  all  pixels  above  it  to  white.  Thresholding 
was  not  used  for  two  reasons.  First,  some  of  the  noise 
would  meet  the  threshold  and  would  be  printed  as  black  as 
the  characters.  This  could  severely  degrade  the  character 
shapes,  causing  the  algorithm  to  fail.  Second,  and  most 
important,  leaving  the  noise  unchanged  was  a  good  test  of 
the  algorithm’s  noise-handling  ability.  It  works  very  well. 

It  may  work  better  without  the  noise;  that  should  be  checked 
in  the  future. 

The  noise  from  the  frozen-frame  effect  can  easily  be 
removed  by  simple  arithmetic-  averaging.  The  noise  between 
any  number  of  scenes  is  uncorrelated  (the  noise  is  assumed  to 
be  white),  yet  the  picture  should  be  highly  correlated.  The 
averaging  process  should  rapidly  reduce  the  noise,  while  not 
reducing  the  scene.  Even  though  a  simple  average  will  remove 
this  noise,  it  again  was  left  it  in  to  see  if  this  algorithm 
would  work  with  the  noise;  it  did. 
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The  source  listing  of  each  program  listed  in  figure  3.2  (page 
16)  and  described  in  chapter  3  is  given  in  this  Appendix. 

The  listing  is  alphabetical,  and  listings  for  any  Fortran 
subroutines  that  the  main  programs  call  are  also  included. 
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Subroutine  CONVERT 


By  Lt.  Simmons 


19  July  81 


This  subroutine  converts  video  pixels  to  lineprinter 
pixels  for  the  Printronix  lineprinter.  This  line- 
printer  requires  that  each  printer  dot  in  a  row  be 
addressed  by  a  bit  of  an  octal  number,  six  dots  per 
octal  number.  A  seventh  bit  must  always  be  one.  The 
eighth  bit  available  in  a  byte  is  zero.  One  peculiar¬ 
ity  is  that  the  bits  addressed  are  reversed  in  the 
byte.  The  right-hand  dot  in  a  row  of  six  dots 
becomes  the  sixth  bit  from  the  right  going  left  in 
the  byte.  The  first  dot  on  the  left  in  the  row  of 
six  dots  becomes  the  first  bit  on  the  right  in  the 
byte.  Note  that  two  pixels  are  processed  in  each 
call  to  CONVERT. 


SUBROUTINE  C0NVERT( IV , ILP) 

DIMENSION  ILP(3)  ;Three  words  per  pixel 

DO  1  11=1,3 

ILP(II)=100K  ;Zero  the  three  words 

IVR=IV. AND. 15  ;Strip  off  right  pixel 

IVL=ISHFT( IV, -4) .AND. 15  ;Strip  off  left  pixel 
IF(IVL.EQ.O)GO  TO  0  ;Left  pixel  is  zero 

GO  TO (5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75)1 VL 
STOP  "Error  in  Subroutine  CONVERT"  ; Should  not  occur 


Process  left  pixel 


ILP(1)=ILP(1) +007 K 
ILP(2)=ILP(2)+007K 
ILP(3)=ILP(3)+007K 
GO  TO  75 

ILP(1)=ILP(1)+007K 
ILP(2)=ILP(2) +00 5 K 
ILP(3)=ILP(3)+001-K- 
GO  TO  75 

ILP(1)=ILP(1) +00 6 K 
ILP(2)=ILP(2)+007K 
ILP(3)=ILP(3) +00 3 K 
GO  TO  75 

ILP(1)=ILP(1) +007 K 
ILP(2)=ILP( 2) +005K 
ILP( 3) =ILP( 3) +002K 
GO  TO  75 

ILP(1 )sILP(1 J+005K 
ILP(2)=ILP(2)+005K 
ILP(3)=ILP(3)+005K 
GO  TO  75 

ILP( 1 )=ILP( 1 ) +00 5 K 
ILP(2)=IIP( 2) +002K 
ILP(3)=ILP(3)+005K 
GO  TO  75 


;Left  pixel  is  0 


;  Jump 

to 

process  right 

pixel 

;Left  pixel  is 

1 

;  Jump 

to 

process  right 

pixel 

;Left  pixel  is 

2 

;  Jump 

to 

process  right 

pixel 

;Left  pixel  is 

3 

;  Jump 

to 

process  right 

pixel 

;Left  pixel  is 

4 

;  Jump 

to 

process  right 

pixel 

;Left  pixel  is 

5 

;  Jump 

to 

process  right 

pixel 

o  o  o 


30 


;Left  pixel  is  6 


ILP(  1)=ILP(1) +00 2K 
ILP(2)=ILP(2)+007K 
ILP(3)=ILP(3)+002K 
GO  TO  75 

35  ILP(1)=ILP(1) +002K 

ILP( 2) =ILP( 2) +005K 
ILP(3)=ILP(3)+002K 
GO  TO  75 

40  ILP(1)=ILP(1)+005K 

ILP(2)=ILP(2)+000K 
ILP(3)=ILP(3)+005K 
GO  TO  75 

45  ILP(1)=ILP(1)+004K 

ILP(2)=ILP(2)+002K 
ILP(3)=ILP(3)+001K 
GO  TO  75 

50  ILP(1)=ILP(1)+002K 

ILP(2)=ILP(2)+000K 
ILP(3)=ILP( 3) +005K 
GO  TO  75 

55  ILP(l)sILPCl) +00 OK 

ILP(2)=ILP(2)+005K 
ILP(3) =ILP(3) +000K 
GO  TO  75 

60  ILP(l)sILPO) +00  4K 

ILP(2)=ILP(2)+OOOK 
ILP(3)=ILP(3)+001K 
GO  TO  75 

65  ILP(1)=ILP(1) +000K 

ILP(2)=ILP(2) +00 2K 
ILP(3)=ILP(3)+OOOK 
GO  TO  75 

70  ILP(1)=ILP(1)+004K 

ILP(2)=ILP(2)+000K 
ILP(3)=ILP(3)+©00K- 

Process  right  pixel 

75  IF(IVR.EQ. 15) RETURN 

IF(IVR.EQ. 14)GO  TO  80 
IFdVR.EQ.  13JG0  TO  81 
IF(IVR.EQ. 12)G0  TO  82 
IFdVR.EQ.  11)GO  TO  83 
IFdVR.EQ.  10)GO  TO  84 
IFdVR.EQ. 9)G0  TO  85 
IFdVR.EQ. 8)GO  TO  86 
IFdVR.EQ. 7)GO  TO  87 
IFdVR.EQ. 6)G0  TO  88 
IFdVR.EQ. 5)G0  TO  89 
IFdVR.EQ.  4)G0  TO  90 
IFdVR.EQ. 3)G0  TO  91 
IFdVR.EQ. 2)GO  TO  92 


;Jurap  to  process  right  pixel 
;Left  pixel  is  7 


;Jump  to  process  right  pixel 
;Left  pixel  is  8 


;Jump  to  process  right  pixel 
;Left  pixel  is  9 


;Jump  to  process  right  pixel 
;Left  pixel  is  A 


;Jump  to  process  right  pixel 
;Left  pixel  is  B 


;Jump  to  process  right  pixel 
;Left  pixel  is  C 


;Jump  to  process  right  pixel 
;Left  pixel  is  D 


;Jump  to  process  right  pixel 
;Left  pixel  is  E 

Corrt-rnue  for  more  processing 


Right 

pixel 

is 

F 

Right 

pixel 

is 

E 

Right 

pixel 

is 

D 

Right 

pixel 

is 

C 

Right 

pixel 

is 

B 

Right 

pixel 

is 

A 

Right 

pixel 

is 

9 

Right 

pixel 

is 

8 

Right 

pixel 

is 

7 

Right 

pixel 

is 

6 

Right 

pixel 

is 

5 

Right 

pixel 

is 

4 

Right 

pixel 

is 

3 

Right 

pixel 

is 

2 
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IF(IVR.EQ. 1)G0  TO  93  jRight  pixel  is  1 

ILP(1)=ILP(1)+070K  ;Right  pixel  is  0 

ILP( 2) =ILP( 2) +070K 

ILP(3)=ILP(3)+070K 

RETURN 

80  ILP(1)=ILP(1)+040K  jRight  pixel  is  E 

ILP(2)=ILP( 2) +000K 

ILP(3)=HP(3)+OOOK 

RETURN 

81  ILP(1)=ILP(1) +000K  jRight  pixel  is  D 

ILP(2)=ILP(2)+020K 

ILP(3) =ILP(3) +000K 
RETURN 

82  ILP(1)=ILP(1)+040K  jRight  pixel  is  C 

ILP(2)=ILP(2)+000K 

ILPC3)=ILP(3)+010K 

RETURN 

83  ILP(1 )=ILP( 1 )+000K  jRight  pixel  is  B 

ILP(2)=ILP(2)+050K 

ILP(3) =ILP(3)+000K 
RETURN 

84  ILP(1)=ILP(1) +020K  *  jRight  pixel  is  A 

ILP(2)=ILP( 2) +00 OK 

ILP(3)=ILP(3)+050K 

RETURN 

85  ILP( 1 )=ILP( 1 )+040K  jRight  pixel  is  9 

ILP(2)=ILP(2)+020K 

ILP(3)=ILP(3)+010K 

RETURN 

86  ILP(1)=ILP(1) +05 OK  jRight  pixel  is  8 

ILP( 2) =ILP( 2) +000K 

ILP(3)=ILP(3)+050K 

RETURN 

87  *  ILP(  1 )  =ILP(  1  )-+D20K  --  jRight  pixel  is  7 

ILP(2)=ILP(2)+050K 

ILP(3)=ILP(3)+020K 

RETURN 

88  ILP(1)=ILP(1)+020K  jRight  pixel  is  6 

ILP(2)=ILP(2)+070K 

ILP(3)=ILP(3)+020K 

RETURN 

89  ILP(1)=ILP(1)+050K  jRight  pixel  is  5 

ILP(2)=ILP(2)+020K 

ILP(3)=ILP(3)+050K 

RETURN 

90  ILP(1)rILP(1) +05 OK  jRight  pixel  is  4 

ILP(2)=ILP( 2) +05 OK 

ILP(3)=ILP(3) +05 OK 
RETURN 

91  ILP(1)=ILP(1) +070 K  jRight  pixel  is  3 

ILP(2)=ILP( 2) +05 OK 
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ILP(3)=ILP(3)+020K 

RETURN 

ILP( 1 )=ILP( 1 )+060K  ;Right  pixel  is  2 

ILP(2) =ILP( 2) +070K 
ILP( 3) =ILP( 3) +030K 
RETURN 

ILP( 1 ) =ILP( 1 ) +070K  ;Right  pixel  is  1 

ILP(2)=ILP(2)+050K 

ILP(3)=ILP(3) +07 OK 

RETURN 

END 


ooo  ooo  ooooooooooooooo 


Program  CREAD 


Version  2 


This  program  will  calculate  the  magnitude  and/or 
phase  of  a  complex  256  by  256  file.  The  output  mag¬ 
nitude  file  is  always  called  "RMAG"  and  the  output 
phase  file  is  always  called  "ANG."  The  input  file 
is  specified  in  the  command  line.  The  command  line 
format  is: 

CREAD/LM/P]  Inputfile 

where  the  global  switch  /M  means  create  RMAG,  and 
the  global  switch  /P  means  create  ANG.  Both  switches 
can  be  be  used.  At  least  one  must  be  given. 

DIMENSION  10(1024) ,101(512) ,102(512) ,IFILE(7) fMS(2) 
COMPLEX  MA(256) 

COMMON  MA,RMAG(256),ANG(256) 

EQUIVALENCEC  MAC  1 ) , I0( 1 ) ) , ( 101 ( 1 ) , RMAGC 1  ) ) ,  (  I02( 1 ) , 

:  ANG(D) 

Accomplish  disk  I/O  tasks. 

CALL  GROUND(I) 

IF( I . EQ . 0 ) OPEN  0 , "COM. CM"  ;0pen  ch.  0  to  "COM. CM" 

IFCI.EQ. 1)0PEN  0 , "FCOM.  CM"  ;0pen  ch.  0  to  "FCOM.CM" 
CALL  COMARGC 0 , I , MS , IER )  ;Read  global  switches 
IFCIER.NE. 1)TYPE"  C0MARG1  error:", IER 
IF(MS( 2) . NE. 0) STOP  "Bad  global  switch." 

IF(MS(1).NE.1. AND. MS ( 1 ) . NE. 8. AND.MSC 1 ) . NE.9)ST0P 
:  "Bad  global  switch." 

IFCMSC 1 ) . EQ. 1 )G0  TO  1  ;Calculate  phase  only 
TYPE"  RMAG  will  be  created." 

DELETE  "RMAG" 

"CALL  CFILWC  "RMAG"  ,-3,512,1) 

IFCI.NE.DTYPE"  RMAG  CFILW  error:", I 

IF( MS( 1 ) . EQ . 8)G0  TO  2  ;Calculate  magnitude  only 

1  TYPE"  ANG  will  be  created." 

DELETE  "ANG" 

CALL  CFILWC "ANG" ,3,512,1) 

IFCI.NE.1)  TYPE"  ANG  CFILW  error:", I 

2  CALL  COMARGC 0 ,IFILE,I,IER) 

IFCIER.NE. 1)TYPE"  C0MARG2  error:", IER 
OPEN  3 , IFILE 

IFCMSC1).NE.1)0PEN  4, "RMAG" 

IFCMSC 1 ). NE. 8) OPEN  5, "ANG" 

Calculate  magnitude  and/or  phase. 

DO  5  1=0,255 

CALL  RDBLKC  3, (1*4), 10, 4, IE) 

IFCIE.NE.1)  TYPE"  RDBLK  error:", IE, I 
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DO  4  J=1 ,256 
RMAG( J ) =CABS( MA( J ) ) 

IF(RMAG(J).EQ.O.O)GO  TO  3 
X=REAL(MA( J ) ) 

Y=AIMAG( MA( J ) ) 

ANG( J ) =ATAN2( Y, X) 

GO  TO  4 

3  ANG( J ) =0 . 0 

4  CONTINUE 

IFCMS(1).NE. 1)CALL  WRBLK( 4, ( 1*2) , 101 , 2, IE) 
IF(MS(  1 ) .  NE.  1 .  AND.  IE.  NE.  1 )  TYPE”  WRBLK4r",IE 
IF(MS(  1 ) .  NE.  8)CALL  WRBLK(  5,  ( 1*2) ,  102, 2, IE) 
IF(MS(1).NE.8.AND.IE.NE. 1)  TYPE  ”  WRBLK5=",IE 

5  CONTINUE 
CALL  RESET 
STOP 

END 
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Program  DISPLAY 
Version  7 


By  Lt.  Simmons 
Fortran  IV 


14  Oct  1981 


This  program  will  convert  video  pixels  to  lineprinter 
pixels,  and  will  put  the  picture  in  a  file  or  on  the 
Printronix  300  lineprinter.  This  program  prints 
either  the  complete  256  by  256  pixel  picture,  or  a 
smaller  picture  that  does  not  contain  the  noise 
created  by  the  video  digitizer  (the  last  five  blocks 
in  the  video  file  are  noise). 

DIMENSION  ILP(3,67) ,IFILE(7) .IRECC64) ,ISAV(3) 

LOGICAL  SHORT 


Solid  line,  space,  and  line  feed/plot-on  characters 


IL=177777K 
NC=40 1 00K 
LF=012K 
LFPC=241 2K 


;Solid  line 
; Space 
;Line  feed 
;Line  feed/plot  on 


Set  up  parameters  for  complete  picture  display. 


SH0RT=. FALSE. 

N1  =66 

N2=256 

N3  =  1 

N4  =  66 

N5  =  67 

N6  =  1 


{Short  picture  test 
;Top  and  bottom  border  length 

{Number  of  lines  displayed 
{Location  of  left  border 
{Location  of  right  border 
{Location  of  line  feed 
{Length  of  lines  displayed 


Open  the  video  file  for  input 

ACCEPT"  What  is  the  name  of  the  input  file?  " 

READ( 1 1 , 17) IFILE( 1 )  {Read  video  file  name 

CALL  0PEN( 1 , IFILE , IER )  {Open  the  video  file 

IF( IER . NE. 1 ) TYPE"  Input  open  error:", IER 


Ask  for  an  output  file,  either  the  lineprinter  or  a 
disk  file,  and  open  the  disk  file  if  necessary. 


ACCEPT"  Do  you  want  a  disk  file  created  (YES/NO)?  " 
READ(11,19)N  {Read  one  ASCII  character 

IF( N, EQ. 19968)GO  TO  2  {File  was  not  selected 
IF( N. NE. 22784)G0  TO  20  {Input  error 

ACCEPT"  A  file  was  selected;  what  should  its  name  be 
(13  char  max)?  " 

READ( 1 1 , 17) IFILE( 1 )  {Read  output  file  name 
CALL  DFILW( IFILE, JER)  {Make  sure  that  the 

IF( JER. NE. 1 .AND. JER. NE. 13)TYPE"  Output  delete  error:" 
,  JER 

CALL  CFILW( IFILE, 2 ,KER)  {file  does  not  exist 


ooo  ooo  ooo  o n o o o  ooo  ooo 


IF(KER. NE. 1 )TYPE"  Output  create  error:", KER 
CALL  OPEN( 12,IFILE,LER)  ;Open  the  output  file 
IF(LER.NE. 1 )TYPE"  Output  file  open  error:", LER 
GO  TO  3 

2  TYPE"  The  picture  will  only  go  to  the  lineprinter. " 
Choose  complete  picture  or  noiseless  picture. 


3 

4 


5 


6 

7 


8 


9 


10 


CONTINUE 

ACCEPT"  Do  you  want  a  complete  picture  (YES/NO)?  " 
READ( 1 1 , 19)N  ;Read  one  ASCII  character 

IF(N.EQ. 19968)GO  TO  22  {Response  was  "NO" 

IF( N. NE. 227  84)G0  TO  21  {Input  error 

TYPE"  A  complete  256  by  256  pixel  picture  was  chosen. 


Put  a  border  at  the  top  of  the  picture. 


IF(SHORT)WRITE( 12, 18) 

DO  7  1=1,3 

IF( SHORT) WRITE  BINARY(12)NC 
IF( SHORT ) WRITE  BINARY(12)NC 
DO  6  J  =  1 , N 1 
WRITE  BINARY ( 1 2) IL 
WRITE  BINARY ( 1 2) LFPC 


{Double  space  first 

{Space  right 
; twice 

{Print  a  line 
{Terminate  the  line 


Each  line  of  the  picture  will  have  a  border  on  each 
end.  A  DO-LOOP  loops  233  (or  256  for  whole  frame) 
times  around  the  next  three  program  parts 


DO  13  JA= 1 , N2 

Put  a  border  down  the  left  hand  side 


DO  8  K=1,3 

'  IF(  SHORT)  ILPCK  ,-1 }  =NC 
IF( SHORT )ILP(K,2)=NC 
ILP(K,N3)=43500K 

Put  a  border  down  the  right 

DO  9  L=1,3 

ILP(L,N4)=40170K 

ILP(L,N5)=LFPC 

Convert  the  video  picture  p 

DO  10  M=1 ,64 

READ  BINARY ( 1 ) IREC(M) 

DO  12  N=N6 , 64 

IWR=IREC(N).AND.377K 
IWL=ISHFT( IREC(N) ,-8) . 
CALL  CONVERT(IWL,ISAV) 


{Put  a  left  border 
;two  spaces  in 
;(~1/4"  space) 

{for  short  picture 

hand  side 

{Insert  border  and 
{line  feed  after 
;the  picture  rows 

ixels  to  $LPT  pixels 


{Read  one  video  line 

{Right  two  pixels 
•377K  {Left  two  pixels 
{Convert  pixels 
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N7=N+1 

IF( SHORT) N7=N 
DO  11  JBsl  ,3 

ILP(JB,N7)=ISHFT(ISAV(JB), 8)  ;Put  in  hi  byte 
CALL  CONVERT( IWR , ISAV)  {Convert  pixels 

DO  12  JC=1 ,3 

ILP( JC,N7)=ILP( JC,N7)+ISAV( JC)  ; Put  low  byte 
DO  13  J  E= 1 , 3 
DO  13  JD=1,N5 

WRITE  BINARY(12)ILP(JE,JD) 

Put  a  border  and  title  at  the  bottom  of  the  picture 
DO  15  JF=1,3 

IF ( SHORT ) WRITE  BINARY(12)NC  {Space  right  twice 

IF ( SHORT) WRITE  BINARY(12)NC  ;for  short  picture 
DO  14  JG=1,N1 

WRITE  BINARY ( 1 2) IL  {Print  a  line 

WRITE  BINARY ( 1 2) LFPC  -.Terminate  the  line 
WRITE  BINARYC 1 2 ) LF  ;End  with  a  line  feed 

WRITE(12,16)  ,-Title  picture 

CALL  RESET  {Close  all  channels 

STOP 

Format  statements. 

FORMAT( *  1 , 15X, 'Signal  Processing  Laboratory,  Air 
Force  Institute  of  Technology,  Wright-Patterson  AFB, 

OH  45433<14>' ) 

F0RMAT(S13)  {Filename  format 

FORMATPO')  {Double  space  for  short  picture 

FORMAT(SI)  {Query  format 

ACCEPT"  Input  error.  Try  again  (YES/NO)  >  " 

GO  TO  1 

ACCEPT"  Input  errors-  Try  again  (YES/NO)  >  " 

GO  TO  4 


Set  up  parameters  for  shortened  picture. 


TYPE"  The  shortened  (noise  removed)  picture  was 


chosen. " 
SHORT=.TRUE. 
N1  =63 
N2=233 
N3=3 
N4=65 
N5=66 
N6  =  4 
GO  TO  5 


{Turn  on  short  test 
{Top  and  bottom  border  length 
{Number  of  lines  displayed 
{Location  of  left  border 
{Location  of  left  border 
{Location  of  right  border 
{Length  of  lines  displayed 


ooo  ooooo 


? 


SUBROUTINE  HISTOGRAM( RMAX, RMIN, IFILE) 


This  subroutine  will  print  out  on  the  terminal 
screen  the  number  of  values  in  25  ranges  across 
the  range  RMIN  to  RMAX. 

DIMENSION  RMAGC256) ,IFILE(7) , NUMBER (25) 

COMMON  IOM(512)  {Integer  array  for  mag.  input 
EQUIVALENCE  ( RMAG( 1 ) , IOM( 1 ) ) 

TYPE"  Subroutine  HISTOGRAM  executing." 

Zero  the  storage  numbers. 

DO  1  1=1,25  ; 25  storage  locations 

1  NUMBERC I) =0 

OPEN  1, IFILE  ;0pen  mag.  file  on  ch.  1 
RINC= ( RMAX- RMIN) *0.04  ; 25  ranges 

DO  2  J  =0 , 255  ; 256  lines 

CALL  RDBLK( 1 , ( J*2 ) ,I0M,2,JER) 

IF( JER. NE. 1 )TYPE"  RDBLK  error:", JER 
DO  2  K=1,256  ; 256  numbers  in  RMAG 

DO  2  L=1,25  ; 25  ranges 

IF((RMAG(K).GE.(RMIN+(L-1)*RINC)).AND. 

( RMAG( K) . LT. ( RMIN+L*RINC) ) ) NUMBER( L) = 

:  NUMBER(L)+1 

IF(NUMBER(L).EQ. 32766 )NUMBER(L) =32765 

2  CONTINUE 

CLOSE  1  {Close  mag.  file 

DO  3  M=1,12  {Two  lines  at  a  time 

0UT1B=RMIN+((M-1)*RINC)  {Bottom  of  1st  range 

0UT1T=RMIN+( M#RINC)  {Top  of  1st  range 

OUT2B=RMIN+( ( M+1 1 )*RINC)  {Bottom  of  2nd  range 

0UT2T=RMIN+((M+12)*RINC)  {Top  of  2nd  range 

N=M+1 2 

3  WRITE(10,4)0UT1B,  OU-TiT,  NUMBER( M) , 0UT2B, 0UT2T, NUMBER( N) 

4  FORMAT ( '  Range  ' , El  1 . 4 , *  to  • , El  1 . 4 , ' : •  ,  15  , '  Range  ', 

:  E 1 1 . 4 , *  to  1 , El  1 .4, » : » ,15) 

WRITE(10,5)0UT2T,RMAX,NUMBER(25)  {Write  25th  range 

5  FORMATC  ',39X, 'Range  ',E11.4,»  to  •  ,  El  1 . 4 ,  •  :  »  ,  15 ) 
RETURN 

END 
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Program  HRV  By  Lt.  Simmons  11  Sep  1981 

Fortran  5 

This  program  energy-normalizes  a  given  magnitude 
array,  then  combines  it  with  a  specified  phase  array 
to  produce  a  256  by  256  array  of  complex  numbers. 

The  command  line  format  is: 

HRV  Magnitudef ile  Phasefile  Outputfile 

INTEGER  CFILE(7) ,PHFILE(7) ,C,PH 
DIMENSION  MFILE(7) , PHASEC 1024) 

REAL  MAG( 1 024 ) , SUMSQ 
COMPLEX  C0UT( 1 024 ) 

COMMON  I0M( 2048 ) ,I0PH( 2048) ,I0C( 4096) 

EQUIVALENCE  (I0M( 1 ) , MAG( 1 ) ) , ( I0PH( 1 ) , PHASE( 1 ) ) , 

:  C  COUT( 1 ) , I0C( 1 ) ) 

M=  1  {Channel  for  magnitude  file 

PH=2  {Channel  for  phase  file 

C=3  {Channel  for  complex  output  file 

Complete  disk  file  tasks. 

CALL  I0F(3,M1 ,MFILE,PHFILE,CFILE, II ,12,13,14) 

OPEN  M, MFILE  {Open  ch.  1  to  mag.  file 

OPEN  PH, PHFILE  {Open  ch.  2  to  phase  file 

DELETE  CFILE 

CALL  CFILW( CFILE, 3 i 1 024 , LER)  {Create  output  file 
IF(LER. EQ. 41 )STOP  "Insufficient  contiguous  blocks: 

:  Output  file" 

IF(LER.NE.1)TYPE"  Output  file  creation  error:", LER 
OPEN  C, CFILE  {Open  output  file  on  ch.  3 

1  FORMAT( S13) ■ 

Read  in  magnitudes  for  energy  normalization. 

SUMSQ=0.0  {Variable  for  storing  sum  of  squares 

Z=65536.  ; 256  squared 

DO  2  1=0,504,8  {Loop  64  times 

CALL  RDBLK(M,I,IOM, 8,NER)  {Read  8  blocks 
IF( NER. N£. 1 )TYPE"  Mag.  RDBLK  error:”, NER 
DO  2  J= 1 , 1024 

SUMSQ=SUMSQ+(MAG( J)/Z)*«2  {Sum  of  squares 

2  CONTINUE 

SUMSQ=SQRT( SUMSQ)  {Square  root  of  sum  of  squares 
REWIND  M  {Rewind  magnitude  file 

Read  in  magnitudes  and  phases  and  calculate  complex 
numbers. 

DO  4  K=0,504,8  {Loop  64  times  again 

CALL  RDBLK (M,K,I0M, 8, I ERR)  {Read  8  mag.  blocks 

IF(IERR.NE.1)TYPE"  2nd  mag.  RDBLK  error:", IERR 
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CALL  RDBLK(PH,K,IOPH, 8, JERR)  ;Read  8  phase  blks. 

IF(JERR.NE.1)TYPE"  Phase  RDBLK  error:", JERR 
DO  3  L=1,102H 

ENM=MAG(L)/SUMSQ  ; Energy-normalize  magnitudes 
X=ENM*COS( PHASE (L) )  ;Real  part 

Y=ENM*SIN( PHASE(L) )  ;Imag  part 

3  COUT(L)=CMPLX(X,Y)  ;Create  complex  number 

CALL  WRBLK(C,(K«2) ,IOC, 16, KERR)  ;Write  complex  #s 
IF(KERR. NE. 1 )TYPE"  WRBLK  error:", KERR 
J»  CONTINUE 

CALL  RESET  ;Close  all  channels 

STOP 

END 
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SUBROUTINE  I0F( N , MAIN , FI ,F2,F3,MS,S1 ,S2,S3) 


Written  by  Lt.  Simmons  10  Sep  1981 

Version  2 

This  FORTRAN  5  subroutine  will  read  from  the  file 
COM. CM  (FCOM.CM  in  the  foreground)  the  program  name, 
any  global  switches,  and  up  to  three  local  file 
names  and  corresponding  local  switches. 

Calling  arguments: 

N  is  the  number  of  local  files  and  switches  to  be 
read  from  (F)COM.CM.  N  must  be  1,  2,  or  3. 

MAIN  is  an  ASCII  array  for  the  main  program  file 
name. 

FI,  F2,  and  F3  are  the  three  ASCII  arrays  to  return 
the  local  file  names. 

MS  is  a  two-word  integer  array  that  holds  any  global 
switches. 

SI,  S2,  and  S3  are  two-word  integer  arrays  that 
hold  the  local  switches  corresponding  to  FI  through 
F3  respectively. 


Dimension  the  arrays. 

DIMENSION  MAIN(7) ,MS(2) 

INTEGER  F1(7) ,F2(7) ,F3(7) ,S1(2),S2(2),S3(2) 

Check  the  bounds  on— N. —  -  ■=- -= ■ 

IF( N. LT. 1 .OR. N. GT. 3) STOP  "N  out  of  bounds  in  IOF." 
Process  the  data  in  (F)COM.CM 

CALL  GROUND(I)  ;Find  out  which  ground  prog,  is  in 
IF( I . EQ . 0)OPEN  0, "COM. CM"  ;0pen  ch.  0  to  COM. CM 
IF(I.EQ. 1)0PEN  0, "FCOM.CM"  ;0pen  ch.  0  to  FCOM.CM 
CALL  C0MARG( 0 , MAIN, MS , IER )  ;  Read  from  (F)COM.CM 
IF(IER.NE. 1)TYPE"  COMARG  error:", IER 
WRITE( 10, 1 )MAIN( 1 )  ;Type  program  name 

1  FORMATd  Program  '  ,  SI  3  ,  ’  running. ' ) 

CALL  COMARG(0,F1 , SI , JER)  ;Read  from  (F)COM.CM 

IF( JER.NE. 1)TYPE"  COMARG  error  (FI):", JER 
IFCN.EQ. 1 )G0  TO  2  ;Test  N 

CALL  C0MARG( 0,F2,S2,KER)  ;Read  from  (F)COM.CM 

IF(KER. NE. 1 )TYPE"  COMARG  error  (F2):",KER 
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M  IF( N. EQ. 2)G0  TO  2  {Test  N 

CALL  COMARG(OfF3,S3,LER)  ;Read  from  (F)COM.CM 

IF(LER. NE. 1 ) TYPE"  COMARG  error  (F3):",LER 
2  CLOSE  0 
RETURN 
END 


oooooo  ooooo 


SUBROUTINE  ITOC(PXLS,OUT) 

This  program  will  convert  the  video  pixels  passed 
in  PXLS  to  complex  numbers  (with  imaginary  parts 
zero),  and  return  them  in  OUT. 

INTEGER  PXLS(1024) 

COMPLEX  OUT( 1024) 

Note  that  PXLS  has  one  dimension  here,  but  is  two- 
dimensional  in  the  calling  program.  A  single  Do-loop 
can  access  all  values  in  a  one-dimensional  array 
much  more  easily  than  a  two-dimensional  array. 

DO  1  1=1,1024 
BIMAG=0.0 

1  OUT(I)=CMPLX(FLOAT(PXLS(I) ) , BIMAG)  ;I  ->  C 

RETURN 
END 
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SUBROUTINE  LOGHISTC RMAXL , RMINL , IFILE) 

This  subroutine  will  print  out  on  the  terminal 
screen  the  number  of  values  in  25  ranges  across 
the  range  RMINL  to  RMAXL,  using  a  logarithmic  scale. 

DIMENSION  RMAG( 256 ),IFILE(7)»  NUMBER ( 25) 

COMMON  I0M(512)  ;Integer  array  for  mag.  input 
EQUIVALENCE  ( RMAG( 1 ) , IOM( 1 ) ) 

TYPE"  Subroutine  LOGHIST  executing." 

Zero  the  storage  numbers. 

DO  1  1=1,25  ; 25  storage  locations 

1  NUMBER( I) =0 

OPEN  1, IFILE  ;0pen  mag.  file  on  ch.  1 
RINC=( RMAXL- RMINL) *0 . 04  ; 25  ranges 

DO  2  J=0 , 255  *,256  lines 

CALL  RDBLK( 1,(J*2),iOM,2,JER) 

IF(JER.NE. 1)TYPE"  RDBLK  error:", JER 
DO  2  K=1,256  ; 256  numbers  in  RMAG 

DO  2  L=1,25  ; 25  ranges 

IF( ( ALOG1 0( RMAG(K) ) . GE. ( RMINL+( L-1 ) *RINC) ) 

:  . AND. ( ALOG1 0( RMAG(K) ) . LT. ( RMINL+L*RINC) ) ) 

:  NUMBER( L) =NUMBER( L) +1 

IF(NUMBER(L).EQ. 32766) NUMBER(L) =32765 

2  CONTINUE 

CLOSE  1  ;Close  mag.  file 

DO  3  M=1,12  ;Two  lines  at  a  time 

OUT1 B=RMINL+( ( M-1)*RINC)  ;Bottom  of  1st  range 

OUT1 T=RMINL+( M*RINC)  ;Top  of  1st  range 

0UT2B=RMINL+( ( M+11)#RINC)  ;Bottom  of  2nd  range 

0UT2T=RMINL+( (M+12)*RINC)  ;Top  of  2nd  range 
N=M+1 2 

3  WRITE ( 10,4)0UT1B,0UT1T,NUMBER(M),0UT2B,0UT2T,NUMBER(N) 

4  F0RMAT( ’  Range  *,E11.4,'  to  * , El  1 . 4 , « : * , 15 , '  Range  ', 

:  Ell. 4  1  to  1  Ell. 4  »:*  15)  _ 

WRITE( 1 0 , 5 )0UT2T  j  RMAXL ) NUMBERC  25 )  ;Write  25th  range 

5  FORMAT ( ’  * , 39X, '  Range  »,E11.4,’  to  » ,E1 1 . 4, * : » ,15) 
RETURN 

END 
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Program  Pick  Version  1  Lt.  Simmons  10  Sep  1981 
This  program  will  display  a  pixel  magnitude  only  if 
it  falls  between  a  user-specified  range.  A  histogram 
(e.g.  program  PIX1)  can  list  possible  pixel  magnitude 
ranges. 

The  command  line  format  is: 

PICK  Inputfile  Outputfile 

INTEGER  FILEK7)  ,FILE0(7)  ,PIX(1024)  ,VIDE0(  256)  ,  SAVE(  4) 
REAL  MAGC1024) 

COMMON  IOM(2048),MIO(2O48) 

EQUIVALENCE  ( I0M( 1 ) , MAG(  1 ) ) 

Perform  disk  I/O  tasks. 

CALL  IOF(2,M1 , FILEI , FILEO , II , M2 , 12 , 13 , 14 ) 

OPEN  1,FILEI  ;0pen  mag.  file  to  ch.  1 

DELETE  FILEO  ;Make  sure  that  FILEO  does  not  exist 
CALL  CFILW( FILEO , 3  » 64 , IER )  ;on  the  disk 

IF(IER.EQ. 41 )ST0P  "Insufficient  contiguous  blocks: 

:  FILEO" 

IF( IER. NE.1) TYPE"  CFILW  error:", IER 
OPEN  2, FILEO 

Request  max.  and  min.  levels  for  pixel  inclusion. 

ACCEPT"  Highest  level  to  include?  " , RMAX 
ACCEPT"  Lowest  level  to  include?  " , RMIN 

Read  all  magnitudes  and  write  magnitudes  that  are  in 
the  specified  range. 

DO  3  1=0,63 

CALL  RDBLK(1,d-*3),IOM,8,JER)  ;Read  8  blocks 
IF( JER.NE. 1 )TYPE"  RDBLK  error:", JER 
DO  1  J=1 , 1024 

PIX(J)=15  ;Make  all  pixels  white 

IF(MAG( J).LE. RMAX. AND. MAG(J).GE. RMIN) PIX(J)= 

:  IFIX(  1 1 . *(MAG(  J )-RMIN)/(  RMAX- RMIN) )  -.Create  pixel 

1  CONTINUE  ;Range  of  pixel  creation  loop 

DO  2  K=0 , 255 

DO  2  L= 1 , 4 

SAVE(L)=PIX( (K*4)+L)  ;Save  4  pixels/loop 
CALL  REPACK ( SAVE,VIDE0(K+1))  ;Repack  4  pixels 

2  CONTINUE  ; Range  of  repacking  loop 

CALL  WRBLK(2,I,VIDE0, 1 ,KER)  ;Write  one  block 
IF(KER.NE.  DTYPE"  WRBLK  error:", KER 

3  CONTINUE  ; Range  of  selection  loop 

CALL  RESET 

STOP 

END 
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Program  PIX1  converts  an  input  256x256  real  number 
file  into  a  VIDEO  file. 

Version  5  2  Sep  1981 

The  command  line  format  is: 

PIX1C/H/L]  Inputfile  Outputfile 

where  the  global  switch  /H  causes  a  histogram  to  be 
printed  on  the  console,  and  adding  the  global  /L 
switch  causes  the  histogram  and  the  clip  level 
requests  to  be  in  base  ten  logarithm  format. 

INTEGER  0UTFILEC7) 

DIMENSION  10(2048) ,101(1024) ,102(256) ,INFILE(7) ,MS( 2) 
COMMON  RMAGO024) 

EQUIVALENCE( RMAG( 1 ) , 10 (  1 ) ) 

Disk  I/O  tasks. 

CALL  I0F(2,M1 , INFILE , OUTFILE , I 1 , 12 , MS , 13 , 14 , 15 , 16) 

IF( MS(  1 ) .  NE.  0 .  AND.  MS ( 1 ) .  NE.  25'6 .  AND.  MS(  1 ) 

: . NE. 272. . OR. MS( 2) . NE. 0)ST0P  "Bad  global  switch.” 

IF( MS( 1 ) . EQ . 272 )TYPE"  Base  ten  logarithm  histogram 
:  option  on.” 

IF( MS( 1 ) . EQ . 256 ) TYPE”  Normal  histogram  option  on." 

OPEN  3, INFILE 

DELETE  OUTFILE 

CALL  CFILW(OUTFILE,3,64,IE) 

IF(IE.NE.  DTYPE"  Output  CFILW:",IE 
OPEN  4, OUTFILE 

Find  min  and  max  values  of  input  file. 

RMIN= 1 . 0E60 
RMAXrO.O 
DO  1  1=0,63 

CALL  RDBLK(3,(8*I),IO,8,IE) 

IF(IE.NE.I)  TYPE"  RDBLK1  error:", IE 
DO  1  J  =  1 , 1024 

RMAX=AMAX1 ( RMAG( J ) , RMAX) 

RMINrAMINI ( RMAG( J ) , RMIN) 

1  CONTINUE 

TYPE"  Min:", RMIN,"  Max:", RMAX 
RMINL=ALOG10( RMIN)  {Logarithm  of  RMIN 

RMAXL=AL0G10(  RMAX)  {Logarithm  of  RMAX 

TYPE"  Min  ( log) : " , RMINL, "  Max  ( log) : " , RMAXL 
REWIND  3  {Rewind  input  file  channel 

Call  histogram  (if  requested). 
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IF(MS(1).EQ. 256) CALL  HISTOGRAM( RMAX , RMIN , INFILE) 
IF(MS(1).EQ.272)CALL  LOGHIST( RMAXL, RMINL, INFILE)  ;Log 

Convert  reals  to  gray  scale  integers  and  pack. 


2 


3 

4 


5 


IF( MS( 1 ).EQ. 256. OR. MS(I).EQ.O) ACCEPT” 
max?  " , RMAX1 

IF( MS( 1 ). EQ . 272 ) ACCEPT”  Log  clip  gray 
RMAX1 

IF( MS( 1 ) . EQ . 272 ) RMAX1 =10 . **RMAX1 
RMAX=AMIN1 ( RMAX, RMAX1 ) 

IF( MS( 1 ) . EQ. 256. OR. MS ( 1 ) . EQ . 0 ) ACCEPT” 
min?  ” , RMIN1 

IF(MS(1).EQ.272)ACCEPT"  Log  clip  gray 


Clip  gray  scale 
scale  max?  ”, 

Clip  gray  scale 
scale  min?  ”, 


RMIN1 

IF(MS(1 ) .EQ.272)RMIN1=10.**RMIN1 
RMIN=AMAX1 ( RMIN , RMIN1 ) 

DO  4  1=0,63 

CALL  RDBLK(3,(I*8),IO,8,IE) 

IF(IE.NE.I)  TYPE”  RDBLK2  error:”, IE 
DO  2  J= 1 , 1024 

RMAG( J) =AMIN1 ( RMAG( J ) , RMAX) 

RMAG( J ) =AMAX1 ( RMAGC J ) , RMIN ) 

I01(J)=IFIX( 15.0*(RMAG( J)-RMIN)/( RMAX-RMIN) ) 
CONTINUE 


KK=0 

DO  3  K  =  1 , 256 
I02(K) =0 
DO  3  J  =  1,4 
KK=KK+1 

I02(K) =ISHFT( 102 (K)  ,4) 
I02(K)=I02(K)+I01(KK) 
CONTINUE 

CALL  WRBLK(4,I,I02,1 ,IE) 

IF(  IE.  NE.  1 )  TYPE-”’  WRBLK  'errdr  :  ” ,  IE 
CONTINUE 

WRITE( 10,5)  OUTFILE(I) 

CALL  RESET 
STOP 

FORMATC  '  ,S13,  'created.  » ) 

END 
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SUBROUTINE  REPACK( N, PIXELS , PXWD) 

Written  by  Lt.  Simmons  Version  2 

This  subroutine  will  repack  four  4-bit  integer  pixels 
into  one  16-bit  word  for  use  by  CHOPS.  Parameter  N 
allows  more  than  one  4-bit  to  1-word  repacking 
operation  in  each  call  to  REPACK. 

INTEGER  PIXELS(4,N) , PXWD(N) 

DO  1  J= 1 , N  ;Loop  N  times 

PXWD( J ) =0 
DO  1  1=1,4 

PXWD( J)=ISHFT(PXWD( J) ,4)  ;Shift  pixel  left  in  word 
1  PXWD( J)=PIXELS(I,J) +PXWDC  J )  ;&  add  next  pixel  on  rt. 
RETURN 
END 
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SUBROUTINE  UNPACK(N,PIXWORD, PIXELS) 


Written  by  Lt.  Simmons  Version  2 

This  subroutine  will  unpack  four  4-bit  integers  from 
a  16-bit  integer  word.  The  pixels  in  a  video  file 
have  to  be  unpacked  if  each  pixel  is  to  be  operated 
on  separately. 

INTEGER  PIXWORD( N) , PIXELS ( 4 , N)  ;Four  pixels  per  word  ' 
DO  1  1  =  1  ,N  ; 1 N *  allows  higher-order 

DO  1  J  =  1,4  {arrays  to  be  passed. 

PIXELS ( (5-J) ,I)=15.AND. PIXWORDC I )  ;Pick  off  rt  pixel 
PIXWORD(I)=ISHFT(PIXWORD(I) ,-4)  ;Shift  word  4  bits 
RETURN  {right  to  pick  off  next  pixel. 

END 
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Program  VIDEO  By  Lt.  Simmons  5  Aug  1981 
Version  2.1 

This  program  sets  up  parameters  for  a  mode  three 
(video)  call  to  the  Cromemco  Z-2D  processor,  via 
the  fortran  subroutine  CHANNEL. 


DIMENSION  IPAR( 2) , IH0LD(7  ) 

INTEGER  CR0ERR,N0VERR,FILE(7) .PDCONT 
LOGICAL  BTEST 

TYPEMNOTICE :  CHOPS  must  be  running ! <1 2>" 

ACCEPT"Input  or  output  ( IN-0/0UT-1 ) ?  ",IDIR 
IF(IDIR. NE.O.AND. IDIR. NE. 1 )G0  TO  8  {Error  check 
1  J=0  {Dummy  parameter 

K=3  {Mode  3  ->  abort 

Send  an  abort  to  the  Cromemco  to  get  its  attention. 


CALL  CHANNELC J, J,K, J, J,«A", J, J, J,IE,IS)  {Abort 


IDCNT=4  ;Set  up 

IPCNT=1  {parameters 

NTSK=3  ; f or  call  to 

IM=2  ; CHANNEL 

ACCEPT"Enter  time  (SEC.):  ",IPAR(1)  ;Display  time 

ACCEPT"What  is  the  name  of  the  data  file  (13  Char 
:  max)?  ” 

READ(11,8)FILE(1)  ; Video  frame  filename 

IF(FILE(1).NE.  10752)G0  TO  3  ,’An  "*”  means  to 

DO  2  1=1,7  ;run  with  the  last 

FILE( I) =IH0LD( I )  ; f ile  used 

IF( IDIR. EQ. 1 )G0  TO  4 

CALL  DELETE(FILE)  {Delete  file  on  input 
TYPE"  Check  monitor” 

CALL  CHANNEL( NTSK, IDIR , IM, IPCNT, IDCNT, FILE ,64,0, 

:  IPAR , IERR , ISYS)  ;Call  to  CHANNEL  — 

IF(IERR.EQ.O.OR. (IERR.EQ. 1 3323 • AND. IDIR . EQ . 0) ) 

:  GO  TO  6  {Jump  if  no  errors 

IF( IERR.EQ. 13323. OR. IERR.EQ. -24 832. OR. IERR.EQ. 

:  -24064. OR. IERR 

:  . EQ.-22528)G0  TO  5  ;Jump  for  specific  error  messages 

IF(BTEST( IERR , 15) )G0  TO  9  {Abnormal  return 

GO  TO  10  {Error  without  abort 

IF( IERR. EQ. -24832)TYPE"<7>ABORT — FILE  DOES  NOT  EXIST" 
IF( IERR.EQ. 13323. OR. IERR.EQ. -24064 )TYPE"<7>AB0RT— 
:FILE  DOES  NOT  CONTAIN  VIDEO” 

IF( IERR. EQ. -22528)TYPE"<7>ABORT — FILE  CANNOT  BE 
:  CREATED" 

TYPE"<12>” 

DO  7  1=1 ,7  ;Save  filename 

IHOLD(I)=FILE(I)  {for  next  loop 

ACCEPT"What  next  ( INPUT-0 , OUTPUT-1 , STOP-2) ?  ”,IDIR 
IF( IDIR. EQ. 0)G0  TO  1 
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IFUDIR.EQ.  1)G0  TO  1 

IF(  IDIR .  EQ .  2)  STOP  ‘'Type  VIDEO  to  rerun.” 
F0RMAT(S13)  {Filename  input  format 

TYPE”<7>Input  error;  try  again.” 

GO  TO  6 


TYPE"<7><12>  ABORT  INITIATED !!!< 1 2>” 
CROERRz 15.AND.IERR 
PDCONT=ISHFT( 240 . AND. IERR , -4) 
NOVERR=ISHFT(-256.AND.IERR,-8) 

TYPE”  CROMEMCO  ERROR  RETURNED CROERR 


There  are 
two  error 
codes  in  the 
variable 
' IERR' 


TYPE”  PARAMETER/DATA  COUNT  RETURNED PDCONT 
CALL  BCLR(N0VERR,7) 

TYPE”  NOVA  ERROR  RETURNED NOVERR  ; Error 
TYPE"  ERROR  CODE  RETURNED IERR  jmessages  are 
TYPE"  DATA  COUNT IDCNT  ; printed  for 

TYPE"  PARAMETER  COUNT IPCNT  ;user  information 

TYPE"  SYSERR  RETURNED ISYS  ;and  correction 

GO  TO  6 


END 


101 


o  o  o  oooonoooooooooooonoooooooooo 


Program  VtoC  **************************************** 

* 

Version  5  * 

* 

Written  by  Lt.  Simmons  22  Sep  1981  * 

* 

*********************************************************** 

This  program  will  read  in  a  video  file  (256  by  256 
pixels)  from  a  disk  and  will  then  write  it  out  to  a 
disk  file  (either  a  contiguous  or  random  file)  as  a 
(256  by  256)  complex  array.  Note  that  the  video  file 
is  only  2 56  by  64,  since  four  pixels  are  packed  in 
each  16-bit  word  in  a  video  file.  These  pixels  have 
to  be  "unpacked”  before  they  can  be  written  out  as  the 
real  parts  of  complex  numbers.  The  first  14  pixels 
in  each  line  can  be  blanked  (made  gray  by  setting  to 
8888  hex.). 


The  command  line  format  is: 

VTOC[/B]  Inputfile  Outputf ile/ t R  or  C] 

where  the  global  switch  /B  invokes  blanking,  and  the 
local  switch  (/R  or  /C)  causes  the  Output  file  to  be 
either  a  random  file  or  a  contiguous  file,  respect¬ 
ively. 

DIMENSION  NAME(7) ,N0UT(7) ,MS(2) ,ISW(2) 

INTEGER  VIDEO( 256 ) , PXLS( 1024)  ; 25 6  words=1024  pixels 
COMPLEX  0UTO024)  ;1024  complex  pixels 

LOGICAL  Z 
COMMON  I0UT(4096) 

EQUIVALENCE  (0UT( 1 ) , I0UT( 1 ) )  ~ 

Z=. FALSE. 

Open  the  video  input  and  output  files. 

CALL  I0F( 2 ,M, NAME , NOUT, II, MS, 12, I SW, 13) 

OPEN  1,NAME  ;0pen  ch.  1  to  video  file 

IF(ISW(1).EQ.8192)G0  TO  2  ;Contiguous  file 

IF( ISW( 2) . NE.  16384)GO  TO  1  ;Error — send  message 

DELETE  NOUT  {Ensure  file  does  not  exist 

CALL  CFILW( NOUT , 2 ,KER)  {Create  random  file 

IF(KER. NE. 1 )TYPE"  Random  file  creation  error:",KER 
GO  TO  3  ; Continue 

1  TYPE"  Bad  local  switch." 

CALL  RESET 
STOP 


102 


ooo  noo  o  o  o  o 


2  TYPE"  Please  wait  while  a  contiguous  file  is  created." 

DELETE  NOUT  {Ensure  file  does  not  exist 

CALL  CFILW( NOUT, 3 ,1024,LER)  {Create  contiguous  file 
IFCLER.EQ.41 )G0  TO  10 

IF(LER. NE. 1 )TYPE"  Contiguous  file  creation  error:" 

:  ,  LER 

C 

3  OPEN  2fN0UT  {Open  the  file  on  ch.  2 

C 

IF(MS( 1 ) . EQ. 0. AND. MS(2).EQ.0)G0  TO  5  ;No  blanking 
IF(MS(1).NE.16384)GO  TO  4  ;Error— send  message 

TYPE"  Gray  blanking  will  occur." 

Z=.TRUE.  ;Blank  option  on 

GO  TO  5 

4  TYPE"  Bad  global  switch." 

CALL  RESET 
STOP 

Process  the  change  from  video  pixels  to  complex 
numbers. 

5  DO  8  1=0,63  ; 6 4  blocks  per  video  file 

11=1*16 

CALL  RDBLK( 1,1, VIDEO, 1,NER)  ;Read  in  1  block 
IF( NER. NE. 1 )TYPE"  RDBLK  error:", NER 
IF((I.GE.59). AND. Z) CALL  ZERO(VIDEO)  {Blank 
IF( . NOT. Z ) GO  TO  7  ;Check  for  blanking 
DO  7  J  =  1 , 193 , 64 

DO  6  K=0,2  {Blanking  steps 

6  VIDEO( J+K) =VIDEO( J+K) . OR .  1 0421  OK  {Blank 

;12  bits 

VIDEO( J+3) =VIDE0( J+3) • OR. 10400K  {Blank  2  bits 

7  CONTINUE 

CALL  UNPACK( 256 , VIDEO , PXLS)  {Unpack  four  pixels 
CALL  IT0C( PXLS , OUT)  -  {Integer  to  complex 

CALL  WRBLK(2,II, IOUT, 16 ,IERR)  {Write  complex  fts 
IF( IERR, NE. 1 )TYPE"  WRBLK  error:", IERR 

8  CONTINUE 
CALL  RESET 
STOP 

Format  statements. 

9  FORMAT(SI)  {Query  format 

Send  error  message  for  insufficient  contiguous  blocks. 

10  TYPE"<7><15>" 

TYPE"  Sorry,  but  there  are  insufficient  contiguous" 
TYPE"  blocks  for  the  output  file  to  be  created." 

TYPE"  Would  you  like  to  stop,  or  run  with  a" 

ACCEPT"  random  file  (STOP/R)?  " 


103 


READOI  ,9)NSR 

IF( NSR. EQ . 21 248) STOP 

IF( NSR. NE. 20992) GO  TO  12 

TYPE"  A  random  file  will  be  created  instead." 

CALL  CFILW(N0UT,2, JERR)  ;Create  random  file 

IF( JERR. NE. 1 )TYPE"  Random  file  creation  error:", JERR 
GO  TO  3 

ACCEPT"  Input  error.  Please  try  again  (STOP/R)  >  " 

GO  TO  11 


O  o  oo  f> 


SUBROUTINE  ZERO(WORD) 

Version  2 

22  Sep  1981 

This  subroutine  gray  blanks  the 

j  j-ariKs  une  word  passed  to  it. 

INTEGER  W0RDC256) 

DO  1  1=1,256 

WORD( I)  =  10421  OK  .c 

RETURN  »Set  w°rd  to  8888  (hex) 

END 
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